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Abstract. 129 Long Period (LP) events, divided into two families of similar events, were 
recorded by the 50 stations deployed on Mount Etna in the second half of June 2008. 
During this period lava was flowing from a lateral fracture after a summit strombolian 
eruption. In order to understand the mechanisms of these events, we perform moment 
tensor inversions. Inversions are initially kept unconstrained to estimate the most likely 
mechanism. Numerical tests show that unconstrained inversion leads to reliable moment 
tensor solutions because of the close proximity of numerous stations to the source po- 
sitions. However, single forces cannot be accurately determined as they are very sensi- 
tive to uncertainties in the velocity model. Constrained inversions for a crack, a pipe or 
an explosion then allow us to accurately determine the structural orientations of the source 
mechanisms. Both numerical tests and LP event inversions emphasise the importance 
of using stations located as close as possible to the source. 

Inversions for both families show mechanisms with a strong volumetric component. 
These events arc most likely generated by cracks striking SW-NE for both families and 
dipping 70° SE (Family 1) and 50° NW (Family 2). For Family 1 events, the crack ge- 
ometry is nearly orthogonal to the dike-like structure along which events are located, 
while for Family 2 the location gave two pipe-like bodies which belong to the same plane 
as the crack mechanism. The orientations of the cracks are consistent with local tecton- 
ics, which shows a SW-NE weakness direction. The LP events appear to be a response 
to the lava fountain occurring on the lO"* of May, 2008 as opposed to the flank lava flow. 
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1. Introduction 

Mt Etna is an active 3330 m high stratovolcano located 
on the East coast of Sicily, Italy. An eruptive period began 
on the 10*'' of May 2008 with a powerful lava fountain in 
the South East Crater, one of the four main summit craters. 
An eruptive fissure opened on the 13"' of May on the East 
flank of the volcano, in the "Valle del Bove" [see e.g., Can- 
nata et al, 2009b]. The flank eruption stopped on July 6*'' 
2009. 

Long Period (LP) events recorded on Mount Etna have a 
period range between 0.75 and 5 s (i.e. 0.2-1.3 Hz). These 
events are commonly defined by a narrow frequency range 
and an emergent onset [see e.g., Chouet, 2003, for a general 
definition of LP events] . In the last few years LP events on 
Mt Etna have been analysed and located in many studies 
[Falsaperla et al, 2002; Saccorotti et al, 2007; Lokmer et al., 
2007, 2008; Patane et al, 2008b; Cannata et al, 2009a]. 
Sources of these events are usually found to be located a 
few hundred meters below the summit craters. LP events 
are repetitive, which suggests a repeating action of the same, 
non-destructive source process. 
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Moment tensor inversion (MTI) has been performed on 
several volcanoes to quantify the source processes of these 
events. Most of these studies [e.g. Ohminato et al., 1998; 
Nakano et al., 2003; Chouet et al., 2003; Kumagai et al., 
2002; Lokmer et al, 2007; Jolly et al, 2010] suggest fluid- 
filled crack mechanisms, often accompanied by single forces. 
Such a mechanism is supported by forward modelling of 
fluid-filled resonator systems with various geometries, such 
as crack or pipe-like conduits [e.g. Chouet, 1985, 1986; Jous- 
set et al., 2004]. This resonance produces slow interface 
waves, also called crack waves, whose dispersive proper- 
ties allow for the generation of low frequency events from 
relatively small sources [Chouet, 1986; Ferrazzini and Aki, 
1987]. The trigger mechanism for this excitation still re- 
mains uncertain, but they are thought to be related to in- 
stabilities in the fluid motion [Ohminato et al, 1998; Rust 
et al, 2008; Neuberg et al, 2006; Cilbert and Lane, 2008]; 
The fluid can be magma [Neuberg et al, 2006], water or 
steam [Cusano et al, 2008; Kumagai et al, 2005], gas [Lok- 
mer et al, 2007], or mixtures of these fluids [Ohminato et al, 
1998]. Considering the hypothesis of a mechanical interac- 
tion between the fluid and the solid phase [Ferrazzini and 
Aki, 1987; Jousset et al, 2003], the strong link between LP 
activity and fluid dynamics implies that the characterisation 
of LP source mechanisms is fundamental to the understand- 
ing of processes in magmatic systems. 

To correctly reconstruct the moment tensor, the topog- 
raphy [Neuberg and Pointer, 2000; Jousset et al, 2004; 
O'Brien and Bean, 2009] and the velocity model [Bean et al, 
2008; Cesca et al, 2008] play important roles. A poor knowl- 
edge of the velocity structure can lead to apparently stable 
but erroneous solutions, as it induces uncertainties in the 
Green's functions. This is particularly true for signals with 
frequencies above 0.2 Hz, whereas the larger wavelengths 
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of the Very Long Period events [Chouet et at, 2003] make 
this approach more stable [Kumagai et al., 2010]. This issue 
can largely be solved by using stations located very close to 
the source positions [Bean et al., 2008; Kumagai et al., 2010]. 

On Mt Etna, the first MTI of LP events was performed 
by Lokmer et al. [2007] complemented by a full investiga- 
tion of the LP properties [Saccorotti et al., 2007; Lokmer 
et al, 2008] and MTI of synthetic data [Bean et al, 2008]. 
These authors suggest that the source mechanism generating 
LP events consist of a subvertical crack striking NNW-SSE, 
with a gas "pulsing" excitation. However their dataset con- 
tains only one station located in the summit area; they sug- 
gested that a larger dataset recorded in the close proximity 
to the sources would help to better constrain the inversion. 
For this reason, a joint Irish (University College Dublin), 
French (Universite de Savoie, Chambery) and Italian (Is- 
tituto Nazionale di Geofisica e Vulcanologia, Catania and 
Pisa) experiment was conducted on Mt Etna in early sum- 
mer during the 2008 eruption. An exceptionally high density 
network of 50 broad-band stations (30 of them were located 
close to the summit), recorded LP events. De Barros et al. 
[2009] located the source positions of 129 selected events be- 
longing to two different families sharing similar waveforms. 
They found shallow source locations in agreement with pre- 
vious studies on Mt Etna [Lokmer et al., 2008; Cannata 
et al., 2009a], but the high resolution locations allowed them 
to determine a temporal migration never observed before. 

In this study we first present the dataset, and then use 
numerical tests to investigate the resolution and robustness 
of constrained [Nakano and Kumagai, 2005; Lokmer et al., 
2007] and unconstrained inversion with the large number of 
stations available here. For real data, unconstrained and 
constrained inversions give similar solutions. The best solu- 
tion is suggested to be a crack mechanism for both families. 
We then comment on the influence of the station distribu- 
tion, noise, mislocation and mismodelling. The interpreta- 
tion of these cracks, striking SW-NE, are finally discussed 
in relation to the volcanic activity. 

2. Data 

Mount Etna is one of the most active volcanoes in the 
world. Since 1995 a period of increased activity began, with 
eruptions occurring approximately every year. Eruptions 
are of two types: intermittent explosive eruptions from the 
summit craters and flank eruptions, with higher effusion 
rates, which originate from lateral fissures. Volcanic activ- 
ity has been monitored since 2003 by a permanent network 
of broadband stations by the INGV Catania observatory. In 
2008, there was a powerful, explosive eruption on the lO**" of 
May. An intense seismic swarm occurred on the IS**" of May, 
and the amplitude of tremors, which are long-duration sig- 
nals related to the eruptive activity with frequency contents 
between 0.5 and 3 Hz on Mt Etna [Patane et al., 2008b], 
strongly increased, just before the opening of an eruptive 
fracture [Napoli et al, 2008; Cannata et al, 2009b]. This 
fissure, striking approximately NW-SE, developed from the 
base of the North East Crater at about 3050 m. a.s.l. and 
rapidly extended downwards to 2500 m. a.s.l. into the slid- 
ing East flank of the volcano in the Valle del Bove. The 
eruption was marked at the beginning by strong Hawaiian 
activity, the effusion rate then decreased to reach a mod- 
erate rate that lasted until the G**" of July 2009. At the 
same time degassing occurred, mainly in the NE crater. LP 
seismic activity was very high before and during the begin- 
ning of the eruption and sharply decreased in June and July 
2008 [Patane et al, 2008a]. The tremor location and the 



geochemistry show the arrival of primitive magma at the 
end of June 2008, without any change in the lava flow rate 
[C'orsaro and Miraglia, 2009; Patane et al., 2008a]. 

A total of 50 stations (including the 16 permanent sta- 
tions) with three-component broadband sensors (30, 40 or 
60 s cut-off period), were installed on Mt Etna between the 
18*'' of June 2008 and the 3'''* of July 2008. In particular, 
30 of them were located at distances shorter than 2 km from 
the summit area (see Figure 1). 

Before analyzing the data, we deconvolve the instrument re- 
sponse from the recorded signals. More than 500 events are 
found using a STA/LTA method on the bandpass (0.2-1.5 
Hz) filtered data. We then classify these events using a cross- 
correlation analysis between all pairs of signals [Saccorotti 
et al., 2007]. We keep the events that give a correlation co- 
efficient greater than 0.9 with all other events on at least 3 
out of the 4 permanent stations close to the summit. We 
obtain 129 events which belong to two families [De Barros 
et al, 2009]. 

The first family (63 events) is only recorded in the first 
two days of the experiment (18*''-19"' of June), while the 
second family (66 events) is distributed over the first four 
days. After the 22"* of June, the amplitudes of the LP 
events decreases by an order of magnitude. In the same 
time period, the amplitudes of the tremors increase. Since 
both LPs and tremors are in the same spectral range, it 
is impossible to extract additional LP events due to poor 
signal-to-noise ratios. Most of the energy of the selected 
events is concentrated between 0.2 and 1.2 Hz, with a peak 
around 0.9 Hz. However, signal spectra show other peaks 
above 2 Hz, and some signals have higher frequency (>10 
Hz) contents. For both families, the waveforms (unfiltered 
and filtered between 0.2 and 1.2 Hz) and the spectral con- 
tents are shown in Figure 2. Although the waveforms are 
quite similar, the spectral peaks are not the same for the two 
families. Family 2 events exhibit a peak frequency slightly 
higher than in the case of Family 1. 

De Barros et al. [2009] located the source of these events 
with a cross-correlation technique. Their analysis was based 
on the similarities between waveforms recorded by the dif- 
ferent stations, which allowed them to measure the time 
delay between each pair of stations. They used a two-step 
technique by first determining an average location from the 
stack data of a given family and then relocating single events 
within that family. The robustness of this approach was 
checked on synthetic data. The source positions are located 
below the summit craters at very shallow depths, between 
and 800m from the summit for the first family and and 400 
m from the summit for the second family. The hypocenter 
positions are clustered into a subvertical, dike shaped struc- 
ture striking NW-SE (Earn. 1) which branches into two 
pipe-like bodies (Fam. 2). The latter structure lies on a 
plane striking SW-NE and dipping 45° NW. Some events 
from the two different families share the same location, thus 
the waveform difference between the two families has to be 
ascribed solely to a different source mechanism. However, 
the similarities of the waveforms indicate a repetitive source 
within each individual family. 

3. Method 

3.1. Moment tensor inversion 
3.1.1. Unconstrained MTI 

We performed a moment tensor inversion in the frequency 
domain as previously used by Auger et al. [2006] and Lok- 
mer et al. [2007]. The n*'' component of the displacement 
field at station s and at the frequency uj, produced by a 
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source located at the position r is denoted «^(r, uj) and can 
be expressed as: 



Fig. 3). Moment tensor components are 1) for a crack: 



= Glp^g{r , uj)Mpq{uj)+G'„p{r , uj)Fp{uj), with n,p,q = x,y,z, 

(1) 

CnpiVjUj) denotes the Green's functions (GF) and 
Gnp^q{r,Lu) their spatial derivatives. We consider single 
forces (-Fp(aj)), and we assume a symmetric moment ten- 
sor (i.e. no rotational effects) with Mpqiuj) — AIqp{uj). 



Equation 1 can then be rewritten in matrix form. The 
data u^(r,aj) are merged in a column vector d and can be 
expressed in a linear form: 



Gm, 



(2) 



where G is the matrix containing the Green's functions and 
their derivatives and m is a column vector of the moment 
tensor components and/or single forces. As we assume that 
the moment tensor is symmetric, only 6 moment tensor com- 
ponents and 3 single forces have to be determined for each 
frequency. The Source Time Function (STF) for each com- 
ponent is then obtained by applying an inverse Fourier trans- 
form. 

The inversion problem is then linear, and equation 2 is 
solved for each frequency by a classical least-square minimi- 
sation. The associated misfit of the waveforms is defined 
by 

(d-Gm)^(d-Gm) 



Mis = 



d^d 



(3) 



As we do not make any a priori assumptions about the mech- 
anisms, the inversion is unconstrained and we can consider 
either the 6 moment components (MT) or the 6 moment 
components and 3 single forces (MT-I-F). In these cases, we 
have 6 or 9 independent parameters to determine for each 
frequency in the frequency band of interest. 

This inversion leads to 6 or 9 independent STFs, which 
need to be decomposed into mechanisms. The Principal 
Component Analysis [Vasco, 1989], which is based on a Sin- 
gular Value Decomposition of the whole set of STFs, allows 
the estimation of a common STF and its contribution to each 
component, that is the scalar moment tensor. However, this 
approach can be quite imprecise, especially if a time shift ex- 
ists between the different components. The scalar moment 
tensor can also be determined using the peak-to-trough am- 
plitude multiplied by the signum function [Chouet et ai, 
2003]. In both cases, the eigenvalues of the 3x3 scalar mo- 
ment tensor give the mechanism and its eigenvectors yield 
the orientation of the principal axes. The solution can then 
be decomposed into the percentage of explosion, CLVD and 
Double Couple [Vavrycuk, 2001]. This last step, the decom- 
position, is however non-unique and unstable. 
3.1.2. Constrained MTI 

Nakano and Kumagai [2005] and Lokmer et al. [2007] con- 
strained the inversion to the particular mechanisms that are 
considered the most likely source mechanisms generating the 
LP events: a Crack (Cr), a Pipe (Pi) and an Explosion 
(Ex). A crack and a pipe are associated with a plane and 
a cylindrical shape, respectively, with the MT components 
{Mxx, Myy, Mzz) = (l,l,3) and (1,2,2) for a Poisson's ra- 
tio of 0.25, respectively. An explosion refers to a purely 
volumetric mechanism with eigenvalues (1,1,1). We use the 
sets of equations given by Nakano and Kumagai [2005] to 
express the Cartesian components of the moment tensor as 
functions of the azimuth angle (j) and dip angle 6 of the sym- 
metry axis (crack normal or longitudinal axis for a pipe, see 
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= A/o (A/^t 4- 2 sin^ 61 cos^ 


Myy 


= Mo(A/^t-|-2sin^6lsin^ 


Mzz 


= Mo lx 1^. + 2 cos^ 6) 
= Mo(sin^6'sin2(/)) 


' MxY 


Mxz 


= Mo (sin 26* cos (/)) 


, Myz 


= Mo (sin 26* sin </i). 



(4) 



2) for a pipe: 



Mxx 
Myy 
Mzz 

MxY 

Mxz 
Myz 



Mo(A//i -I- cos^ 6)008^ (f) + sin^ 0) 
Mo (A/^i -I- cos^ 6 sin^ (f> + cos^ 0) 
Mo {\/^i + sm^ 9) 
-l/2Mo(sin^6'sin 20) 
-l/2Mo(sin 29 costf)) 
-l/2Mo(sin 29 smtf)), 



(5) 



and 3) for an explosion 



Mxx = Myy = Mzz 
MxY = Mxz = Myz 



Mo 




(6) 



Ado denotes the seismic moment, A and fj, are the Lame's pa- 
rameters. X, Y and Z refer to the East, North and vertical 
upward direction, respectively. Herein, azimuth (f) and dip 9 
are defined using the convention given in figure 3, i.e. is 
measured between and 360° anticlockwise from East and 
9 is defined between and 90° from the upward direction. 
It is worth stressing that the convention we use differs from 
the one used by Nakano and Kumagai [2005]. 

For each frequency, equation 2 can now be rewritten, as: 
d = GMmM(A/Ai,6',<;/>)Mo-hGFmF, (7) 

where Gm are the GF derivatives associated with the mo- 
ment tensor components and niM is the vector containing 
the moment tensor components as defined in equations 4, 5 
and 6, and Mo is the Source Time Function (STF) in the 
frequency domain. Gf and niF are the Green's functions 
and the source properties associated with the single forces, 
respectively. The last term of this equation refers to the 
inversion for single forces, and can either be included or 
omitted. If omitted, i.e. forces are not considered, inversion 
is denoted by Cr, Pi and Ex. As the vector niM is indepen- 
dent of the frequency, for given values of A//i, 9, <f) and uj, 
the inversion procedure reduces to an inversion for a single 
parameter, Mo{(jj). We search for the most likely solution by 
performing a grid search over the 9-(j> domain and making 
an inversion for the STF for each 9-(j} pair. If single forces 
are considered (inversion denoted Cr-|-F, Pi-fF and Ex+F), 
the moment tensor components are still determined using 
a grid search over the 9-(l> domain, while keeping the inver- 
sion for the forces unconstrained, as we do not assume any 
relationship between the single forces and the moment com- 
ponents. For given values of A/^, 9, cj> ^nd uj, the number 
of unknowns to be determined is 3 for the single forces plus 
one for the moment components. 



3.2. Green's function (GF) computations 

3.2.1. 3D Elastic algorithm 

The Green's functions are computed using the 3D Elas- 
tic Lattice algorithm of O'Brien and Bean [2004], based on 
a discrete particle scheme. The model is three-dimensional 
and includes topography of Mt Etna. It is centered on the 
volcano summit and has an area of 19.2x16x7 km with a 40 
m grid size. Absorbing boundaries (4.8 km wide) are applied 
at the bottom and the edges of the model in order to pre- 
vent refiections from the model boundaries. The absorbing 
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boundaries could be smaller in this study, but are designed 
to also work efficiently in the VLP frequency range, i.e. up 
to 0.02 Hz. As the topography strongly distorts waveforms 
[Neuberg and Pointer, 2000; Jousset et al., 2004; O'Brien 
and Bean, 2009], a free surface based on the Digital Eleva- 
tion Model (DEM) of Mt Etna is required. The source func- 
tion used for the GF computation is a Gaussian pulse with a 
7.5 Hz cut-off frequency giving a fiat spectrum below 1.5 Hz, 
as the signals are filtered in the 0.2-1.2 Hz range. Below 1.5 
Hz, the source signature is similar to the spectral signature 
of a Dirac delta function, which should theoretically be used 
to compute Green's functions. We have 25 grid points per 
minimum wavelength up to 1.5 Hz to prevent numerical dis- 
persion for modelling including topography [Ohminato and 
Chouet, 1997]. 

3.2.2. Velocity model and source locations 

Bean et al. [2008] show that the moment tensor is very 
sensitive to incorrect velocity models, and particularly to 
shallow, low velocity structures; they also show that the ef- 
fect of an erroneous velocity model is stronger for stations 
further from the source. In this study, since 30 stations are 
located in the source near-field (less than one wavelength 
away) and we do not have any information on the shallow 
velocity properties, we choose to use a homogeneous model. 
Velocities for P— and 5*— waves are asummed to be 2000 m 
s~^ and 1175 m s~^, respectively. These velocities are simi- 
lar to the results of the recent tomographic study of Mount 
Etna [Monteiller et al., 2009], which is homogeneous in the 
shallow part of the volcano, and to those determined in the 
location process of the LP events considered in this study 
[De Barros et al., 2009]. For stations close to the source, as 
shown by Lokmer and Bean [2010], the near-field effect has 
a strong influence on waveforms, but it is fully taken into 
account in our simulations. Although attenuation is strong 
in volcanic structures [Jousset et al, 2004], it is not as im- 
portant as scattering and topographic effects [O'Brien and 
Bean, 2009]. As the propagation distance is less than one 
wavelength and attenuation is also unknown in the shallow- 
est part of the volcano, it is not considered here. 

Chouet et al. [2003] show the importance of a correct 
source location. However, as shown by Lokmer et al. [2007], 
in the presence of a poorly resolved shallow velocity model, 
there is a trade-off between the source location and source 
mechanism, i.e. coupled inversion can lead to an erroneous 
solution. Moreover, the GF calculation for multiple sources 
with such a large number of receivers is computationally 
expensive, for both direct and reciprocal approaches [Eis- 
ner and Clayton, 2001], as GFs have to be computed from 
each source location or each receiver position. Hence we 
use the source location from De Barros et al. [2009]. The 
events do not share exactly the same source location, but 
we want to avoid the computation of the GFs for multi- 
ple sources. Moreover, the location can be slightly wrong 
because of the trade-off between the location and the veloc- 
ity model. We assume an average source position for both 
families described in section 2 and we use these positions 
to compute the GFs required for the inversion of all LP 
events. In the following sections, we investigate the errors 
introduced by this assumption in synthetic tests and on real 
data. The average source positions have UTM coordinates 
of (499.4, 4178.76, 2.84) for Family 1 and (499.5, 4178.45, 
3.0) km for Family 2, i.e. 490 and 330 m below summit level 
(see Figure 1). 



4. Inversion of Synthetic Data 
4.1. Numerical tests 

In the previous section, we introduced two strong assump- 
tions which are i) a homogeneous model, and ii) an average 
source location for all events. To assess the sensitivity of our 



inversion to uncertainties in the velocity structure, source 
mislocation and noise, we perform inversion of noisy nu- 
merical data computed with a velocity model and a source 
location different from those used in the Green's functions 
calculation. We also intend to assess: 1) if a constrained 
inversion gives more reliable results than an unconstrained 
inversion, and 2) if single forces have to be considered in the 
inversion. In the following sections, when dealing with nu- 
merical results, we use the notation X, Y and Z to refer to 
the East, North and vertical upward directions, respectively. 

The velocity model used to compute the synthetic data is 
an artificial model consisting of a gradient with a Vp increase 
from 1600 m s^^ to 2500 m s~^ from the surface to 500 m 
below the topography. The source location is misplaced by 
90 m in the horizontal plane and 120 m vertically downward 
compared to the position where the GFs are computed. The 
source function is a Ricker wavelet, with 1 Hz central fre- 
quency. This analysis is carried out using synthetic data 
computed for two mechanisms: 1) Vertical crack, called Cx 
with eigenvalues Mxx = 3Mo, Myy = Mo, Mzz = Mo 
(see eq. 4) and the crack-normal oriented along A"— axis 
((^=0°and 61=90°) and amplitude Mo = 3.10^^ N.m; 2) 
Same crack Cx with an added single force F^^ with compo- 
nents {Fx,Fy,Fz) = {9,9,9^2)10^ N, which corresponds 
to an orientation of = 45° and ^=45° . Random noise is 
generated in the same spectral range as the waveforms. The 
amplitude of the noise added to the synthetics is chosen to 
be 25 % of the maximum amplitude at etsm station, in order 
to achieve a noise level rougly similar to the one present in 
the real data. 



4.2. Results 

For the data computed using a crack Cx without single 
forces, the mechanism and its orientation are reconstructed 
well by unconstrained inversion with, and without, single 
forces (see Figure 4 and Table 1). The reconstructed mo- 
ment however has a higher amplitude than the true solution. 
This is an expected result as i) the GF source is deeper (120 
m) than the data source location and ii) the low velocity 
near the surface amplifies the signal relative to the GFs, 
calculated in a higher velocity. We then invert the synthetic 
data for a Crack, Pipe and Volumetric constrained mecha- 
nism (see Table 1). The minimum residual is obtained for 
the correct mechanism, i.e. for the crack constrained inver- 
sion. Its orientation and the STF are correctly recovered. 
As the solutions with and without noise are similar (not 
shown in figure), MTI does not appear to be very sensitive 
to the noise in this case. 

The wrong velocity model and the mislocation produce 
spurious forces, which contribute to about 50 % in the re- 
construction of the waveforms. This can be due to i) the 
distortion of the waves propagating in the inhomogeneous 
velocity model and ii) difference in arrival times for all sta- 
tions when waveforms are computed with mislocation and 
wrong velocity models. We observe that moment tensors do 
not seem to be strongly affected by these factors for such 
a near- field deployment. The errors in the modelling seem 
to be accommodated by the single forces, i.e. the inversion 
leads to spurious single forces to reconstruct the data when 
the modelling is not perfect. Moreover, the spurious forces 
are very similar for both inversions presented in Figure 4, 
as they are representative of the model errors. In conclu- 
sion, i) constrained, or not, the inversion does not help to 
determine if single forces are real or due to mismodelling in 
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the GF computations, and ii) unconstrained inversions are 
reliable for reconstructing the MT part of the source mech- 
anisms for such a station network. 

The same inversions are carried out for the second dataset 
where the true source is a crack Cx and a strong single 
force _p45 (see Fig. 5 and table 1). The inversion for mo- 
ment only (MT) gives an incorrect mechanism, similar to 
a pipe (1,1.8,2.2). The moment tensor components recon- 
structed by the unconstrained inversion (MT-fF) are very 
close to the true solution. The orientation of the main axis 
is however close to the true solution for both cases (less than 
15° error in dip and azimuth). For crack constrained inver- 
sions (Cr and Cr-|-F), moment is not reliable when forces 
are not considered as the reconstructed crack appears to 
be horizontal instead of vertical, while the inversion con- 
strained for a crack and single forces (Cr-|-F) gives very good 
results. Among the different geometries (Pi, Cr and E), the 
misfit minimum is not necessarily obtained for the correct 
mechanism (see table 1). As shown by Lokmer and Bean 
[2010], radiation patterns for Crack and Pipe mechanisms 
are very similar, and can appear the same if the wavefield 
sampling is not sufficiently dense. When including single 
forces, the misfit minimum leads to a crack mechanism, but 
the misfit differences are still very small. Similarly to the 
first case, the forces found in these inversion tests are not 
properly reconstructed as they include spurious forces due 
to velocity mismodelling and mislocation. The difference 
in misfit values between inversion with and without single 
forces is larger than for case 1, but it is not large enough to 
assure that forces are real. In conclusion, forces cannot be 
accurately reconstructed, but inversion with single forces is 
required if single forces exist. 

In order to make our results more general, we consider 
additional synthetic data computed for different crack ge- 
ometry and single forces. The geometry of the crack called 
Cz, is inspired by the results of Lokmer et al. [2007] and is de- 
fined hy 4> = Zb° ,e = 72° and amplitude Mo = 3 10^^ N.m. 
We also consider a vertical force, designed by Fz with com- 
ponents (0,0,6 10^) N, as several authors [e.g. Chouet et al., 
2003] found a single force with this orientation. Data are 
then computed with the same conditions as in the previous 
tests, which include mislocation, mismodelling and noise. 
Results are summarized in table 1, for the different combi- 
nation of the two cracks and the two forces. In agreement 
with the previous examples, the MTI is more accurate when 
single forces are included, as the solution obtained with the 
MT-I-F inversion is closer to the true mechanism than the 
solution reconstructed by the MT only inversion. The ori- 
entation of the major axis is however quite stable whatever 
the unconstrained inversion we used, the errors on the az- 
imuths and dips are less than 20°. Moreover, the Cr-|-F 
inversion leads to a very good solution in all cases, while the 
Cr solution is sometimes totally wrong. When single forces 
are not included, the misfit cannot be used to discriminate 
between the different geometries in constrained inversions, 
as the solution for pipe constrained inversion can have a 
smaller misfit than the one for a crack. 



4.3. Conclusion and strategy 

In the case of mismodelling and mislocation, moment ten- 
sor inversions do not allow single forces to be properly recon- 
structed. However, numerical tests show that the moment is 
more reliable if the inversion is carried out considering free 
single forces. Consequently, herein we allow single forces 
in the inversion to compensate for the errors coming from 
the velocity model and the source location. These single 



forces are however not considered for the interpretation of 
the mechanism. A similar conclusion has been reached by 
Sileny [2009], who shows that, for earthquakes with double- 
couple mechanisms with a small non-shear component, and 
in the presence of mislocation and mismodelling, solutions 
are more stable when considering the 6 MT components than 
for a constrained double-couple inversion. Therefore, it can 
be sometimes better to have more unknowns in an inver- 
sion process in order to accommodate the errors. However, 
synthetic tests, like the ones presented here, are always nec- 
essary to choose an inversion strategy (constrained MTI?, 
single forces?), as results will strongly depend on station 
density and topography. 

With their stations configurations. Bean et al. [2008] 
showed that the strong sensitivity of the MTI to the shallow 
velocity model prevents the recovery of the correct uncon- 
strained solution and can lead to totally wrong orientation of 
the mechanism. However, in the study presented here, as the 
mismodelling and mislocation effects seem to be balanced by 
using stations very close to the source, it is possible to cor- 
rectly reconstruct the mechanism using an unconstrained 
inversion. The interpretation of such a solution can be deli- 
cate, as the extraction of a single STF and the decomposition 
of the moment tensor components can be delicate and im- 
precise. To help alleviate this problem, Chouet et al. [2005] 
and Kumagai et al. [2005] first determine a rough approxi- 
mation of the mechanism using an unconstrained inversion, 
and then refine the source characteristics by constraining 
the inversion. We choose to use a similar approach, which 
proceeds in two steps: 1) we invert for an unconstrained 
solution in order to determine the most reliable mechanism 
type (Crack, Pipe, Explosion), and 2) we use results of step 
1 to perform a constrained inversion. The second step allows 
us to confirm and check consistency of this solution and to 
refine the mechanism and its characteristics. In both steps 
we consider forces and leave them unconstrained in order to 
accommodate the errors. 



5. Moment Tensor Inversion of Mt Etna 
Data 

We invert data from both families to determine the mo- 
ment tensor using 16 stations (see section 6 for a justification 
of the number of stations) with the best azimuthal distri- 
bution and signal-to-noise ratio. Since some stations were 
not available at the beginning of the experiment, the set of 
stations is different for each family. Individual events are 
contaminated by noise and do not share exactly the same 
source position [De Barros et al, 2009]. The GFs are how- 
ever computed for fixed positions (see section 3). To check 
this assumption (see discussion in section 6.2), we carried 
out inversion for 44 and 39 events for Family 1 and 2, re- 
spectively. The mean STF is obtained by averaging all the 
reconstructed STFs and the standard deviation gives us er- 
rors associated with noise and mislocation. To do that, 
although the LPs have very similar amplitude, we normalise 
the STFs by the maximum of Mxx in order to have com- 
parable solutions for the different LPs. Errors for Family 1 
(fig. 6) are larger than for Family 2 (Fig. 7), but the cal- 
culated STFs do not show any strong variations for either 
family. As expected for LP multiplets, the source process 
is highly repetitive. Constrained inversion is then carried 
out using a single event in order to show that inversion of a 
single event is also stable. 



5.1. Family 1 

For Family 1, the STF reconstructed by unconstrained in- 
version (16 stations) with and without single forces are very 
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similar (fig. 6a and b). Tlie misfit value (see table 2) is how- 
ever considerably lower when forces are considered. Forces, 
whether physical or an artefact, do not change the moment 
tensor solution in this case. Waveform matches between 
data and reconstructed waveforms are shown in Figure 8 for 
the 16 stations used for the unconstrained inversion with 
single forces for an individual event. Fits are very good for 
most of the stations very close to the source position. They 
disimprove for stations with lower amplitude signals due to 
the lower signal-to-noise ratio, and because the inversion 
gives more weight to the signals with the largest amplitude. 

The unconstrained inversions (MT and MT-I-F) lead to 
very similar mechanism and orientation. They can be in- 
terpreted as a crack (e.g., 1,1.1,2.3 for MT-I-F inversion). 
This is confirmed by the eigenvectors of the moment tensor 
solution, which are shown in figure 9. In the second step, we 
invert for a crack solution, with and without single forces, 
and search for the azimuth and dip (Fig. 6c and d) . Whereas 
STFs are slightly different for the four inversion results, they 
show a very similar amplitude and orientation of the crack 
mechanism. These results are in close agreement with the 
numerical tests shown in Figure 4. We can therefore assume 
that the single forces are probably not real or are too weak 
to be reconstructed. The moment tensor components show 
a crack whose normal is oriented with azimuth <j)=-40° and 
dip 6=70° . 



5.2. Family 2 

For Family 2, we also used 16 stations to perform MTI. 
The moment tensor solutions of the inversions with and 
without single forces are very different, both for the mech- 
anism and for the STFs (see Figure 7). The misfit differ- 
ence between these two inversions is very large (see Table 
2). The waveforms cannot be properly explained without 
forces, though they are very well reconstructed when forces 
are considered (see Figure 10) . By analogy with the numer- 
ical tests shown in Figure 5, we are more confident with the 
solution reconstructed with single forces. However, the time 
shifts which exist among the different moment components 
and single forces do not allow for an easy interpretation of 
the mechanisms. If we use a Principal Component Analysis 
(PC A) for the moment part [Vasco, 1989], the first princi- 
pal component shows eigenvalues of (1,1.1,1.6), with nearly 
80% isotropic component. The deviatoric part shows an 
axis pointing in the {(f)=110° , 8=50°) direction. As the non- 
diagonal components of the tensor are shifted compared to 
the isotropic part and have a weaker amplitude, their effects 
are therefore underestimated by the PCA approach. If we 
use only the maxima of the STFs, we find a source mecha- 
nism of (1,1.4,2.2) with the same orientation for the major 
axis. To have a clearer idea of the mechanism, we plot the 
eigenvectors of the MT-I-F solutions (Fig. 9). The eigenvec- 
tors are more spread out than for Family 1 but they however 
show a clear pattern for the mechanism, with a longer axis 
and two smaller axes with comparable length. This clearly 
indicates a mechanism similar to that of a crack. 

To reconfirm the mechanism, we perform constrained in- 
versions for a crack, a pipe and an explosion (see table 2). 
Smaller misfit values are obtained for the crack mechanism. 
Similarly to the unconstrained inversion, this solutions lead 
to a major axis in the 0=110° and ^=50° direction. The pipe 
and the explosive constrained inversions show slightly higher 
misfits. We choose to discard these mechanisms as they do 
not show a major axis with an orientation consistent with 
the one obtained by the unconstrained inversion. This case 
is very similar to the inversion of synthetic data computed 
for a crack with a single force (Fig. 5). Part of these forces 
may be real but a significant portion is a consequence of 
model error. As shown by synthetic tests, it is impossible to 
accurately reconstruct the forces. For this reason, the single 
forces will not be quantitatively described and discussed. 



6. Discussion 

6.1. Station distribution and density 

The major part of Family 1 events were recorded by only 
16 stations. In order to have comparable results we chose 
to perform the inversion for both families using 16 stations. 
However, Family 2 events and part of the Family 1 signals 
are recorded by 30 stations. We investigate the influence 
of the number of stations, azimuth repartition and propa- 
gation distance, by inverting a single event from Family 2 
with different numbers of stations. For a number of sta- 
tions n between 4 and 28, we randomly choose up to 120 
sets of stations, and perform an unconstrained MT inver- 
sion. We also compute the average propagation distance 
and the azimuthal coverage. The latter is estimated by di- 
viding the azimuthal space into n angular sectors with an 
angle of 360/n° , and calculating the number of portions in 
which at least one station is included. This value is then 
normalized by n to have an estimation of the azimuth cov- 
erage quality as a percentage. The solution obtained using 
the 30 stations is taken as a reference, as this is the best 
solution we can expect from this dataset. We then look at 
the RMS error between the 6 moment components obtained 
using n stations and the whole set of stations. Results are 
summarised in figure 11. We see that no matter how many 
stations are used in the inversion, it is still possible to have 
a solution close to the reference. However, the likelihoods 
vary significantly and at least 10 stations are required to 
be sure that the solution is not totally divergent from the 
reference. Propagation distance is clearly the crucial point, 
as the solutions diverge from the reference when the average 
distance increases. Sensitivity to the azimuth seems to be 
weaker; the decrease of the misfit with azimuthal coverage 
is less systematic. It is however never very bad, at least half 
of the azimuthal sectors are usually covered. Finally, the 
solution shown in figure 6 and computed with 16 stations is 
very similar to the reference solution. 

In conclusion, using between 10 and 30 stations we find 
that the number of stations does not change the solution 
very significantly as long as: 1) a broad azimuthal distribu- 
tion is achieved and 2) some stations closest to the sources 
are used. This is mainly because the stations closest to the 
source strongly constrain the solution. Using less than 8 
stations, MTI usually leads to a different and certainly er- 
roneous solution. This result, however, cannot be directly 
generalised to any other station distributions and any other 
volcano, but it should be tested using synthetic modelling 
for each individual case. Moreover, we do not use the same 
set of 16 stations for Family 1 and Family 2. We also verify 
that the difference between the solutions obtained for the 
two families is not produced by the station distribution dif- 
ference. Inversion of Family 2 events leads to a similar result 
when using either set of stations. 

Lokmer et al. [2007] show that the stations closest to the 
source help to constrain the source time function, while the 
others can be used to determine the mechanism. Here, the 
STF does not resemble the waveform recorded at the clos- 
est station (summit station etsm). Signals from this station 
show a complex waveform (see Fig. 10) probably due to 
local site effects and strong topographic effects. However, 
the other stations with small offset from the source display 
a signal very similar to the STF. In general, for stations 
close to the source, LP waveforms are not strongly distorted 
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by propagation effects, which stabilises the STF reconstruc- 
tion. In conclusion, MTI requires stations in close proximity 
to the source (near summit stations in our case) to be accu- 
rate, but not necessarily a dense network of stations (i.e. a 
minimum of 10 near summit stations in our case). 

6.2. Reliability of the solutions 

The moment tensor inversion solutions for both families 
suggests a crack mechanism. Figure 12 graphically summa- 
rizes the orientation of the solution for the moment compo- 
nents, with the orientation of the main location structures 
found by De Barros et al. [2009]. We are confident of the 
crack solution found for Family 1 as it appears stable for all 
of the inversion tests. For Family 2 the unconstrained solu- 
tion shows a high volumetric component with a mechanism 
between a crack and an explosion. The constrained inver- 
sions and the eigenvector plot confirm the crack solution, 
therefore we have confidence in this solution. 

The principal moment component for both families are 
about (1,1,2). To obtain these values for a crack, we need 
a Poisson's ratio of i/ = 1/3, which implies A = 2/i. This 
high ratio is classically related to the high temperature in 
volcanic rocks [Chouet et al., 2003]. Fractured and uncon- 
solidated media can also present a high Poisson's ratio [e.g. 
Bourbie et al., 1986], which is most likely the case in the near 
subsurface of volcanoes. Another possible explanation of 
the difference with the theoretical crack mechanism (1,1,3) 
is that the latter is computed for an idealised point source 
with no realistic boundary conditions between the cracks 
perimeter and the surrounding medium. We also note the 
presence of strong single forces, especially for Family 2. If 
some of these forces are real they could be related to mass 
movement or drag forces [Takei and Kumazawa, 1994]. It 
also could mean that the mechanism is more complex than 
can be solved by the MTI used here as: 1) the source can 
comprise of several time-delayed mechanisms, 2) it can have 
a spatial extent with complex geometry, and 3) it involves 
rotational effects, such as torque. In these cases, mech- 
anisms cannot be described by a second-order symmetric 
moment tensor and require a higher order moment tensor 
or a different approach to be accurately determined. Nu- 
merical experiments have yet to be conducted to investigate 
how to recover complex sources such as pressure dipoles, 
torque effects, etc. 

Small time shifts between the different components of the 
solution and strong single forces can be partly explained by 
the source mislocation. For the events used in this study 
the common source locations are associated with standard 
errors on the location of 150 m and 65 m for Fam. 1 and 
2., respectively [De Barros et al., 2009]. These errors are 
smaller than the one used for the synthetic tests (175m), 
which does not prevent the reconstruction of accurate MT 
solutions. Furthermore, the standard deviations of the so- 
lutions (Fig. 6 and 7) are very small; this shows that the 
solution is not strongly affected by the choice of the common 
source locations. Mismodelling of the velocity structure is 
another reason for the strong single forces and time shifts. 
Lateral velocity variations have not been considered, either 
in the synthetic test or in the GF modelling. They can how- 
ever be strong, notably because of the deconsolidation of 
the East fiank of the volcano, which is sliding. The velocity 
contrast is certainly stronger in the shallowest part of the 
volcano, which can lead to greater errors for the shallowest 
family (Fam. 2). Finally, a third explanation for the errors 
in the solution can be the complexity of the source; how- 
ever, being limited by the accuracy of the velocity model, it 
is impossible to unambiguously solve the issue of the source 



complexity itself. 



6.3. Source mechanisms 

The orientations of the crack normals are {(f)— -40° , 
0=70°) and (0=110°, 61=50°) (see Figure 12), which cor- 
respond to cracks striking in the SW-NE and SWS-NEN 
directions. Uncertainty of the inversion is ±10°. Cracks 
of both families are roughly orthogonal but their strikes 
show a similar orientation (between N40° E and N70°E). 
The orientation of these cracks are different from the crack 
found by Lokmer et al. [2007] {cj>=35° and 9=72°). How- 
ever, those events were recorded during the 2004 eruption 
and showed different waveforms and spectral characteristics. 
As expected, the LP seismicity on Mt Etna is not constant, 
both for the seismicity rate and for the source properties 
[Patane et al, 2008b]. 

For Family 1, the MT solution (crack) is roughly orthogo- 
nal to the source area determined by De Barros et al. [2009] , 
which is a subvertical dike-like structure striking NW-SE 
(see Fig. 12). For Family 2, the crack solution lies on the 
plane (with normal oriented by (/)=120°and ^=45°) contain- 
ing the two pipe-like bodies. The LP events are not neces- 
sarily produced by the structure in which they are located. 
As expected, the inversion of all the individual events (see 
Fig. 6 and 7) shows that the mechanism is highly repetitive. 

For both families, the source time function is very short 
(i.e. less than 4 s), which suggests a pulsing rather than 
an oscillating mechanism. The STFs can be seen as an im- 
pulsive function filtered in the frequency band of interest 
(0.2-1.2Hz). Amplitudes of the seismic moment are about 
43 and 25 10® N.m for Family 1 and 2, respectively. Volu- 
metric change AV^ can be estimated from Mq = fiAV. From 
the velocity of the medium we compute a rough approxima- 
tion of the dynamic shear modulus: fi = 2.9 GPa, which 
we assume to be equivalent to the static one. The volume 
changes are 15 and 9 m"^ respectively. These volumes are 
smaller than the one found by Lokmer et al. [2007] , but they 
are in agreement with the lower amplitude of the signals and 
the shallower source positions. These volumes correspond to 
a normal displacement of 1 mm for a 100 meter sided square 
crack and 10 cm for a 10 m one. 



6.4. Relationship to the eruption 

De Barros et al. [2009] and this study show that 1) two 
families of LP events are spread out along structures located 
between 800 m and the surface, 2) their source mechanisms 
are related to cracks, with orientations which are not nec- 
essarily the same as the structures along which the events 
are located, 3) for both families, crack strikes are roughly 
similar while dips are orthogonal, and 4) signals and source 
mechanisms are similar within each family. 

Following the lava fountain of the lO"* of May 2008, an 
eruptive fissure opened on the eastern flank of the volcano. 
During the experiment time, lava flowed below the LP source 
locations at a distance of 1.5 km. Because only gases were 
being released from the summit craters, it is possible that 
no magma had risen to the upper part of Mount Etna. This 
suggests that LP events are maybe not directly related to 
magma movements. Moreover, gases may be the most likely 
fluids present in the main conduits and in the fractures sur- 
rounding them. 

The lava fountain at the beginning of the eruption was 
associated with high fluid pressure that can destabilize the 
edifice by opening fractures in the upper part of the vol- 
cano. After this event, the flank lava flow and the summit 
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degassing certainly drained the cone producing a decrease 
of the pressure. The LP events may be Unked to this de- 
compressive phase. The decrease of pressure can lead to 
instabilities of those fractures as the volcano settled under 
its own weight. Patane et al. [2008b] analysed a family of 
LP events (called Family 2 in Patane et al. [2008b]) occur- 
ring only after the lava fountains of 2007, which have similar 
characteristics to the events studied here. They lasted for 
approximately one month after the lava fountain and were 
interpreted as the response to the volcano deflation. This is 
also confirmed by Falsaperla et al. [2002] who linked the LP 
activities to the collapses of the crater floor. In this case, 
the fluids involved can be gas or steam. Gases contained in 
the cracks are suddenly expelled to the main conduits. This 
can produce LP events with mechanisms similar to hydraulic 
transients [Ferrick et al, 1982] or hydrodynamic instabili- 
ties of nonlaminar fiows [Rust et al., 2008]. This hypothesis 
can be linked to the laboratory studies performed by Ben- 
son et al. [2008], who show that the decompression phase in 
rock samples can generate LP events in complex-shape frac- 
tures belonging to the damage zone of the main conduits. 
Similarly to our study, the events are located in the main 
conduits but are not directly generated inside them. 

Cracks for both families are striking SW-NE. The orien- 
tation of the cracks and of the location structures are con- 
sistent with the tectonic setting, which generate faults in 
the NW-SE and SW-NE direction [Bonnacorso and Davis, 
2004]. However, as they are in the shallowest part of the 
volcano, they are more likely due to gravity effects. In 
particular, the East flank of the volcano is collapsing and 
successive eruptions strongly destabilized this area. This 
generates weaknesses oriented SW-NE, which are coherent 
with the azimuth of the cracks. The cessation of the LP 
events after the 22"'' of June suggests that the upper part 
of the volcano reached an equilibrium, where pressure and 
stress return to a static state. The decompression phase, 
following the lava fountain of the 10*'' of May, 2008, lasts 
about 40 days, which is in agreement with the conclusion of 
Patane et al. [2008b]. 



7. Conclusion 

Two families of LP events, comprising 63 and 66 events 
respectively, are selected from the first four days of a seis- 
mic experiment on Mount Etna (18/06/2008-03/07/2008). 
50 stations, including 30 stations located in close proximity 
to the summit, were used for this study. Moment tensor in- 
version of numerical data shows that, for this deployment, it 
is more reliable to use forces in the inversion to correctly de- 
scribe the moment. However, the forces cannot be correctly 
reconstructed as they strongly reflect the errors coming from 
the velocity mismodelling and source mislocation. In gen- 
eral, as MTI appears strongly sensitive to station distribu- 
tion, numerical tests are therefore required before every MTI 
study. Stations close to the source positions are required to 
correctly invert Long Period events. 

We perform moment tensor inversions in two steps. First, 
we determine the type of mechanism involved using uncon- 
strained inversion. We then constrain the inversions to this 
particular mechanism to confirm the solution and refine its 
characteristics. Inversions of the events of the two fam- 
ilies show mechanisms with high volumetric components, 
most likely generated by cracks. For both families, cracks 
are striking in the SW-NE direction, while their dips are 
approximately orthogonal to each other. The crack orienta- 
tions are thus different from the location structures obtained 



by De Barros et al. [2009]. This can suggest that the LP 
events are generated by the fractures which belong to the 
damaged zone around the main conduits of the volcano. We 
hypothesize that these events are related to the decompres- 
sion phase following the lava fountain of the 13"* of May, 
and not to the lava fiow from the fiank eruption. 

MTI reveals strong forces, especially for Family 2, but we 
are not able to determine if they are real or due to artefacts 
in the moment tensor inversion. These forces, as well as the 
time shifts observed between the moment tensors compo- 
nents, can be due to uncertainties in the velocity model or 
to complex sources (e.g., dual sources, torque) that cannot 
be accurately reconstructed by second-order MTI as used 
here. To solve this problem, a better knowledge of the ve- 
locity model is needed. To be able to unambiguously explain 
both moment and forces, a more general approach must take 
into account: i) extended sources, ii) multiple sources, and 
iii) rotational effects. 

Acknowledgments. L. De Barros and G.S. OBrien were 
part funded by the department of Communications, Energy and 
Natural Resources (Ireland) under the National Geosciences pro- 
gramme 2007 2013. Financial assistance for field work from Uni- 
versity College Dublin and the EU VOLUME project, and compu- 
tational facilities from SFI/HEA ICHEC are acknowledged. We 
wish to thank M. MoUhoff, A. Braiden, J. Grangcon, P. Bas- 
cou, M. La Rocca, D. Galluzzo and S. Rapisarda for assistance in 
the field experiment. The relevant comments made by the asso- 
ciate editor and the reviewers Philip M. Benson, Jan Sileny and 
Philippe Joussct greatly contributed to improve the manuscript. 



References 

Auger, E., L. D'Auria, M. Martini, B. Chouet, and P. Daw- 
son, Real-time monitoring and massive inversion of source pa- 
rameters of very long period seismic signals: An application 
to Stromboli Volcano, Italy, Geophys. Res. Lett., 33{i), doi: 
10.1029/2005GL024703, 2006. 

Bean, C, I. Lokmer, and G. O'Brien, Influence of near-surface 
volcanic structure on long-period seismic signals and on mo- 
ment tensor inversions: Simulated examples from Mount Etna, 
J. Geophys. Res., 113, doi:10.1029/2007JB005468, 2008. 

Benson, P. M., S. Vinciguerra, P. G. Meredith, and R. P. 
Young, Laboratory Simulation of Volcano Scismicity, Science, 
522(5899), 249-252, 2008. 

Bonnacorso, A., and P. Davis, Modeling of ground deformation 
associated with recent lateral eruptions: mechanics of magma 
ascent and intermediate storage at Mt Etna, in Mt Etna: Vol- 
cano Laboratory, pp. 293-306, AGU, 2004. 

Bourbie, T., O. Coussy, and B. Zinszner, Acoustics of porous me- 
dia, Institut Frangais du petrole, Paris, France, 1986. 

Cannata, A., M. Hellwcg, G. D. Grazia, S. Ford, S. Alparone, 
S. Gresta, P. Montalto, and D. Patane, Long period and very 
long period events at Mt. Etna volcano: Characteristics, vari- 
ability and causality, and implications for their sources, J. Vol- 
canol. Geotherm. Res., iS7(3-4), 227-249, 2009a. 

Cannata, A., P. Montalto, E. Privitera, G. Russo, and 
S. Gresta, Tracking eruptive phenomena by infrasound: May 
13, 2008 eruption at Mt. Etna, Geophys. Res. Lett., 36, doi: 
10.1029/2008GL036738, 2009b. 

Cesca, S., J. Battaglia, T. Dahm, E. Tcssmer, S. Heimann, and 
P. Okubo, Effects of topography and crustal heterogeneities on 
the source estimation of LP event at Kilauea volcano, Geophy. 
J. Int., 172(3), 1219-1236, 2008. 

Chouet, B., Excitation of a buried magmatic pipe: a seismic 
source model for volcanic tremor, J. Geophys. Res., 90, 1881- 
1893, 1985. 

Chouet, B., Dynamics of a fluid-driven crack in three dimen- 
sions by the flnite difference method, J. Geophys. Res., 91, 
doi:10.1029/JB091iB14pl3967, 1986. 

Chouet, B., Volcano seismology, Pure appl. geophys., 160, 739- 
768, 2003. 



DE BARROS ET AL.: LP EVENT SOURCE MECHANISM AT MT ETNA 



X - 9 



Chouet, B., P. Dawson, T. Ohminato, M. Martini, G. Saccorotti, 

F. Giudiccpietro, G. De Luca, G. Milana, and R. Scarpa, 
Source mechanisms of explosions at Stromboli Volcano, Italy, 
determined from moment-tensor inversions of vory-long-period 
data, J. Geophys. B,es., 108, 2019-2042, 2003. 

Chouet, B., P. Dawson, and A. Arciniega-Ceballos, Source mech- 
anism of vulcanian degassing at Popocatepetl volcano, Mexico, 
determined from waveform inversions of very long period sig- 
nals, J. Geophys. Res., 110{B7), doi:10.1029/2004JB003524, 
2005. 

Corsaro, R., and L. Miraglia, Dynamics of magma in the plumb- 
ing system of Mt. Etna volcano, Sicily, Italy: A contribution 
from pctrologic data of volcanics erupted from 2007 to 2009, 
AGU, Fall Meeting 2009, abstract #V51C-1690, 2009. 

Cusano, P., S. Petrosino, and G. Saccorotti, Hydrothermal ori- 
gin for sustained Long-Period (LP) activity at Campi Flegrei 
Volcanic Complex, Italy, J. Volcanol. Geotherm. Res., 177(A), 
1035-1044, 2008. 

De Barros, L., C. Bean, I. Lokmer, G. Saccorotti, L. Zuccarello, 

G. O'Brien, J. -P. Metaxian, and D. Patane, Source geometry 
from exceptionally high resolution long period event observa- 
tions at Mt Etna during the 2008 eruption, Geophys. Res. 
Lett, 36, doi:10.1029/2009GL041273, 2009. 

Eisner, L., and R. W. Clayton, A Reciprocity Method for 
Multiple-Source Simulations, Bull. Seism. Soc. Am., 91(3), 
553-560, 2001. 

Falsaperla, S., E. Privitera, B. Chouet, and P. Dawson, Anal- 
ysis of long-period events recorded at Mount Etna (Italy) in 
1992, and their relationship to eruptive activity, J. Volcanol. 
Geotherm. Res., 114(3-4), 419-440, 2002. 

Ferrazzini, V., and K. Aki, Slow waves trapped in uid-Ued innitc- 
crack: implication for volcanic tremor., J. Geophys. Res., 92, 
doi:10.1029/JB092iB09p09215, 1987. 

Ferrick, M. G., A. Qamar, and W. F. St. Lawrence, Source mech- 
anism of volcanic tremor, J. Geophys. Res., 87, 8675-8683, 
1982. 

Gilbert, J. S., and S. J. Lane, The consequences of fluid motion 
in volcanic conduits. Geological Society, London, Special Pub- 
lications, 307(1), 1-10, 2008. 

Jolly, A. D., S. Sherburn, P. Joussot, and G. Kilgour, Eruption 
source processes derived from seismic and acoustic observa- 
tions of the 25 September 2007 Ruapehu eruption-North is- 
land. New Zealand, J. Volcanol. Geotherm. Res., 191(1-2), 
33-45, 2010. 

Jousset, P., J. Neubcrg, and S. Sturton, Modelling the time- 
dependent frequency content of low-frequency volcanic earth- 
quakes, J. Volcanol. Geotherm. Res., 128, 201-223, 2003. 

Jousset, P., J. Neuberg, and A. Jolly, Modelling low-frequency 
earthquakes in a viscoelastic medium with topography, Geo- 
phy. J. Int., 159, 776-802, 2004. 

Kumagai, H., B. Chouet, and M. Nakano, Waveform inversion of 
oscillatory signatures in long-period events beneath volcanoes, 
J. Geophys. Res., 107(B11), 2301-2312, 2002. 

Kumagai, H., B. A. Chouet, and P. B. Dawson, Source process 
of a long-period event at Kilauea volcano, Hawaii, Geophy. J. 
Int., 161(1), 243-254, 2005. 

Kumagai, H., et al.. Broadband seismic monitoring of active vol- 
canoes using deterministic and stochastic approaches, J. Geo- 
phys. Res., 115, B08,303, doi:doi:10.1029/2009JB006889, in 
press, 2010. 

Lokmer, I., and C. J. Bean, Properties of the near-field term 
and its effect on polarisation analysis and source locations of 
long-period (LP) and very-long-period (VLP) seismic events 
at volcanoes, J. Volcanol. Geotherm. Res., In Press, doi: 
10.1016/j.jvolgeores. 2010.02. 008, 2010. 

Lokmer, I., C. J. Bean, G. Saccorotti, and D. Patano, Moment- 
tensor inversion of LP events recorded on Etna in 2004 using 
constraints obtained from wave simulation tests, Geophys. Res. 
Lett, 34, doi:10.1029/2007GL031902, 2007. 

Lokmer, I., G. Saccorotti, B. D. Licto, and C. J. Bean, Tempo- 
ral evolution of long-period seismicity at Etna Volcano, Italy, 
and its relationship with the 2004-2005 eruption. Earth Planet. 
Sci. Lett, 34, 205-220, 2008. 

Monteiller, V., J.-L. Got, D. Patanc, G. Barberi, and O. Cocina, 
Double-difference tomography at Mt Etna volcano: Prelimi- 
nary results, in The VOLUME Project. VOLcanoes: Under- 
standing subsurface mass moveMEnt, edited by C. J. Bean, 
A. K. Braiden, I. Lokmer, F. Martini, and G. S. O'Brien, 2009. 



Nakano, M., and H. Kumagai, Waveform inversion of volcano- 
seismic signals assuming possible source geometries, Geophys. 
Res. Lett, 32, doi:10.1029/2005GL022666, 2005. 

Nakano, M., H. Kumagai, and B. A. Chouet, Source mechanism 
of long-period events at Kusatsu-Shirane Volcano, Japan, in- 
ferred from waveform inversion of the effective excitation func- 
tions, J. Volcanol. Geotherm. Res., 122(3-4), 149-164, 2003. 

Napoli, R., G. Currenti, C. Del Negro, F. Greco, and D. Scan- 
dura, Volcanomagnetic evidence of the magmatic intrusion on 
13th May 2008 Etna eruption, Geophys. Res. Lett, 35, doi: 
10.1029/2008GL035350, 2008. 

Neuberg, J., and T. Pointer, Effects of volcano-topography on 
seismic broadband waveforms, Geophy. J. Int., 143, 239-248, 
2000. 

Neuberg, J. W., H. Tuffen, L. Collier, D. Green, T. PoweU, and 

D. Dingwell, The trigger mechanism of low-frequency earth- 
quakes on montserrat, J. Volcanol. Geotherm. Res., 153(1-2), 
37-50, 2006. 

O'Brien, G., and C. Bean, A 3D discrete numerical elastic lat- 
tice method for seismic wave propagation in heterogeneous 
media with topography, Geophys. Res. Lett, 31(1A), doi: 
10.1029/2004GL020069, 2004. 

O'Brien, G. S., and C. J. Bean, Volcano topography, structure 
and intrinsic attenuation: Their relative influences on a simu- 
lated 3d visco-elastic waveflold, J. Volcanol. Geotherm. Res., 
183(1-2), 122-136, 2009. 

Ohminato, T., and B. A. Chouet, A free surface boundary con- 
dition for including 3d topography in flnite-difference method. 
Bull. Seism. Soc. Am., 87, 494-515, 1997. 

Ohminato, T., B. A. Chouet, P. Dawson, and S. Kedar, Wave- 
form inversion of very long period impulsive signals associated 
with magmatic injection beneath Kilauea Volcano, Hawaii, J. 
Geophys. Res., 103, doi:10.1029/98JB01122, 1998. 

Patane, D., A. Cannata, G. D. Grazia, E. Giampiccolo, P. Mon- 
talto, E. Privitera, and Gruppo di Analisi di Sismologia, 
Quadro di sintesi e aggiornamento al 16 luglio 2008 suUo stato 
di attivit sismica dell'etna, INGV report, 1st. Naz. di Geofis. 
e Vulcanol., Rome., 2008a. 

Patane, D., G. Di Grazia, A. Cannata, P. Montalto, and 

E. Boschi, The shallow magma pathway geometry at 
Mt. Etna volcano, Geochem. Geophys. Geosyst, 9, doi: 
10.1029/2008GC002131, 2008b. 

Rust, A. C, N. J. Balmforth, and S. Mandre, The feasibiUty of 
generating low-frequency volcano seismicity by flow through a 
deformable channel, Geological Society, London, Special Pub- 
lications, 307(1), 45-56, 2008. 

Saccorotti, G., I. Lokmer, C. Bean, G. D. Grazia, and D. Patane, 
Analysis of sustained long-period activity at Etna Volcano, 
Italy, J. Volcanol. Geotherm. Res., 160, 340-354, 2007. 

Takei, Y., and M. Kumazawa, Why have the single force and 
torque been excluded from seismic source models?, Geophy. J. 
Int., 118(1), 20-30, 1994. 

Vasco, D., Deriving source-time functions using principal compo- 
nent analysis. Bull. Seism. Soc. Am., 79, 711-730, 1989. 

Vavrycuk, V., Inversion for parameters of tensile earthquakes, J. 
Geophys. Res., 106(B8), 16,339-16,355, 2001. 

Sileny, J., Resolution of non-double-couplc mechanisms: Simula- 
tion of hypocenter mislocation and velocity structure mismod- 
ehng, Bull. Seism. Soc. Am., 99(4), 2265-2272, 2009. 



C. J. Bean, School of Geological Sciences, University College 
Dublin, Belfield, Dublin 4, Ireland. 

L. De Barros, School of Geological Sciences, University College 
Dublin, Belfield, Dublin 4, Ireland, (louis.debarros@ucd.ie) 

I. Lokmer, School of Geological Sciences, University College 
Dublin, Belfield, Dublin 4, Ireland. 

J.-P. Metaxian, LGIT, Univorsite de Savoie-IRD-CNRS, 73376 
Chambery, France. 

G. S. O'Brien, School of Geological Sciences, University Col- 
lege Dublin, Belfield, Dublin 4, Ireland. 

D. Patanc, INGV-Catania, Piazza Roma 2, 95125 Catania, 
Italia. 

G. Saccorotti, INGV-Pisa, Via della Faggiola 32, 56126 Pisa, 
Italia. 

L. Zuccarello, INGV-Catania, Piazza Roma 2, 95125 Catania, 
Italia. 



X - 10 



DE BARROS ET AL.: LP EVENT SOURCE MECHANISM AT MT ETNA 



Table 1. Misfit (Mis) in percentage and moment eigenval- 
ues (MTE) for different constrained (Cr, Pi, Ex denote crack, 
pipe and explosion constrained inversion, respectively) or un- 
constrained (MT) inversion with and without single forces (F). 
The true mechanism involved vertical (Cx) and inclined (C^) 
cracks, as well as vertical (Fz) and inclined (-F45) single forces. 



True Mech. 




Unconstrained 


Crack 


Pipe 


Explosion 


MT 


F 




MT 


MT-fF 


Cr 


Cr-fF 


Pi 


Pi-f-F 


Ex 


Ex-fF 


Cx 


/ 


Mis 
MTE 


29 
1:1.2:2.4 


24 
1:1.1:2.8 


46 
1:1:3 


28 
1:1:3 


52 
1:2:2 


47 
1:2:2 


69 
1:1:1 


63 
1:1:1 


Cx 


F45 


Mis 
MTE 


47 
1:1.8:2.2 


27 
1:1.4:3.4 


63 
1:1:3 


37 
1:1:3 


56 
1:2:2 


39 
1:2:2 


74 
1:1:1 


54 
1:1:1 


Cx 


Fz 


Mis 
MTE 


30 
1:1.6:2.1 


21 
1:1.4:3.2 


61 
1:1:3 
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1:1:3 
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1:2:2 
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1:2:2 
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Cl 
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1:2:2 


25 
1:2:2 
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1:1:1 


31 
1:1:1 
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1:1.1:2.4 
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Table 2. Misfit value (Mis), moment tensor eigenvalues (MTE) for both families and different inversion. 



Family 


Inv 


Mis (%) 


MTE 




MT+F 


27 


1:1.1:2.2 


Fl 


MT 


61 


1:1.5:3.2 




Cr+F 


38 


1:1:3 




Cr 
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1:1:3 




MT+F 


21 


1:1.1:1.6 




MT 


67 


1:3.1:3.8 


F2 


E+F 
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1:1:1 


Pi+F 


40 


1:2:2 




Cr+F 


38 


1:1:3 




Cr 


84 


1:1:3 
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Figure 1. Broadband station positions on Mt Etna. Left: Mt Etna location (top left panel) and all 
stations available on Mt Etna between the 18*'' of June and the 3'^'' of July 2008 (bottom left panel). 
Contour interval is 250 m. Middle: Summit area of the volcano with stations located within 2 km from 
the summit. Contour interval is 100 m. Stars indicate the average LP source location for both families 
found by De Barros et al. [2009]. Right: Same locations (stars) in a North-South cross-section at an 
UTM longitude of 499.5 km (marked by the dashed line in the middle panel). 




Figure 2. Data of an event from a) Family 1 and b) Family 2, recorded at et81 station, vertical com- 
ponent. Top panels: waveforms (raw data and filtered data between 0.2 and 1.5 Hz) and lower panels: 
spectral content. 
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Figure 3. Coordinate system used to define crack and pipe orientation. 
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Figure 4. Source Time Function (black tliin lines) reconstructed using a) unconstrained inversion for 
moment only (MT), b) unconstrained inversion for moment and single forces (MT+F), c) constrained 
inversion for crack only (Cr), d) constrained inversion for crack and single forces (Cr+F). Noisy synthetic 
data are computed for the vertical crack Cx without single forces, whose source mechanism is shown by 
the thick grey lines. In the figure, amplitude is 10^^ N.m for the moment and lO' N for the single forces. 
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Figure 5. Source Time Function (black thin lines) reconstructed using a) unconstrained inversion for 
moment only (MT), b) unconstrained inversion for moment and single forces (MT-I-F), c) constrained 
inversion for crack only (Cr), d) constrained inversion for crack and single forces (Cr-|-F). Noisy synthetic 
data are computed for the vertical crack Cx and inclined single force F^^, whose source mechanism is 
shown by the thick grey lines. In the figure, amplitude is lO^'^ N.m for the moment and 10^ N for the 
single forces. 
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Figure 6. MTI results for Family 1: Mean solution and standard deviation for the Source Time Func- 
tion obtained by the unconstrained inversion of 44 events for a) moment only (MT) and b) moment and 
single forces (MT+F); Source Time Function obtained by constrained inversion of a single event for c) 
crack constrained inversion (Cr); d) crack constrained inversion with single forces (Cr+F). Amplitude is 
10^ N.m for the moment and 10*^ N for the forces. 
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Figure 7. MTI results for Family 2; Mean solution and standard deviation for the Source Time Func- 
tion obtained by the unconstrained inversion of 39 events for a) moment only (MT) and b) moment and 
single forces (MT+F); Constrained inversion of a single event for c) crack constrained inversion (Cr); d) 
crack constrained inversion with single forces (Cr+F). Amplitude is 10'' N.m for the moment and 10^ N 
for the forces. 
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Figure 8. Waveform (displacement) fit between the data (daslied lines) and the synthetic seismograms 
(continuous lines) for an individual event of Family 1. Inversion is unconstrained and includes single 
forces (MT+F). The misfit value is 27 % (see Table 2). 





Figure 9. Eigenvectors obtained for the MT+F inversion solutions for a) Family 1 and b) Family 2. 
Corresponding solutions are given in figure 6b for Fam. 1 and 7b for Fam. 2. Eigenvectors are sam- 
pled every 0.04 s when the amplitude of one of the components is greater than 50% of the maximum 
amplitude. 
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Figure 10. Waveform (displacement) fit between the data (dashed lines) and the synthetic seismograms 
(continuous lines) for an individual event of Family 2. Inversion is unconstrained and includes single 
forces (MT+F). The misfit value is 21 % (see tab. 2). 
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Figure 11. RMS errors of the source time function versus the number of stations, for a single event 
from Family 2. The reference is the solution obtained when using the full set of 30 stations. The line 
show the mean. The stars refer to the inversion results shown in figure 6 and 7. Colorscale corresponds 
to a) azimuthal repartition in % and b) average source-receiver distances in metre. 
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Figure 12. Source mechanisms obtained for a) Family 1 and b) Family 2 events, for the crack con- 
strained inversion (Cr+F). The left panels are map views, and the right panels are 3D views. The circular 
areas represent the cracks. The solid lines represent the normalized eigenvectors, the longest lines are the 
crack normals. The light grey squares show the location structures obtained by De Barros et al, [2009], 
i.e. a) a sub- vertical dike striking SWS-NEN for Family 1 and b) a 45° inclined plane striking SW-NE 
containing the two pipe-like bodies of Family 2. The dashed lines are the normal vectors of these two 
structures. The azimuths of the normal vectors are given on the map view. 



